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Abstract 

Let A = {«!, 02, ... , as} and {bi, 62, • • • , h} with s + r = n, the many to many 
point matching with demands and capacities matches each point G A to 
at least a, and at most a'i points in B, and each point bj G i? to at least Pj 
and at most P'j points in A for all 1 < i < s and 1 < j < t. In this paper, 
we present an O(n^) time and 0{n) space algorithm for this problem. 
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1. Introduction 

A matching between two sets defines a relationship between their ele- 
ments. The matching is used in various fields such as computational biology 
1|, pattern recognition computer vision [sj, music information retrieval 
4|, and computational music theory js]. A many-to-many matching between 
A and B assigns each point in A to one or more points in B, and vise versa. 

Let A and B be two sets with \A\ + \B\ = n, Eiter and Mannila |6i] pro- 
posed an 0{n^) algorithm for the minimum many-to-many matching problem 
between A and B by reducing the problem to the minimum-weight perfect 
matching problem in a bipartite graph. 

The minimum many-to-many matching with demands and capacities, here 
called MMDC matching, is a matching in which each point G A is matched 
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to at least and at most a'i points in 5, and each point hj G -B is matched 
to at least fij and at most _B' ^ points in A, such that sum of the matching costs 
is minimized. Schrijver [7| solved the MMDC matching problem in strongly 
polynomial time. In this paper, we present a new algorithm that computes 
an MMDC matching between A and B in O(n^) time using 0{n) space. In 
section [21 we review the basic Hungarian algorithm and some preliminary 
definitions. In section [31 we present our new algorithm. 

2. Preliminaries 

Given an undirected bipartite graph G = [AUB.E], K maximum match- 
ing M is a matching that for any other matching M', we have Weight{M') < 
Weight{M). A path with the edges alternating between M and E — M is 
called an alternating path. Each vertex v that is incident to one edge in 
M is called a matched vertex; otherwise it is a free vertex. An alternating 
path that its both endpoints are free is called an augmenting path. Note 
that if the M edges of an augmenting path is replaced with the E — M 
ones, its size increases by 1. Let V = A U B, a. vertex labeling function 
I : V ^ R assigns a label to each vertex v & V . A vertex labeling that in 
which /(a) + l{h) > Weight{a, b) for all a G A and 6 G -B is called a feasible 
labeling. The equality graph of a feasible labeling I is a graph G = (V, Ei) 
such that El = {{a, b)\l{a) + l{b) = Weight{a, b)}. The neighbors of a vertex 
u E V is defined as Ni{u) = {v\{v,u) G Ei}. Consider a set of the vertices 
S CV, the neighbors of S is Ni{S) = [j^^^^Ni{u). 

Lemma 1. Consider a feasible labeling I of an undirected bipartite graph 
G = {AU B,E) and S C A with T = Ni{S) ^ B, let 

ai = min {/(flj) + l{bj) — Weight{ai,bj)}. 

If the labels of the vertices of G is updated such that: 

{l{v) — ai if V E S 
llv) + ai if VET 
l{v) Otherwise 

then, V is also a feasible labeling. 

Proof. Note that I is a feasible labeling, so we have l{a) + l{b) > 
Weight{a, b) for each edge (a, b) of E. After the update four cases arise: 



• a G S" and b G T. In this case 

l'{a) + I'ip) = l{a) -ai + l{h) + ai = l{a) + l{h) > Weight{a, b). 

• a ^ S and b ^ T. We have 

l'{a) + I'ib) = l{a) + l{b) > Weight{a, b). 

• a ^ S and b & T. We see that 

l'{a) + l'{b) = l{a) + lib) + ai> 1(a) + l{b) > Weight{a, b). 

• a E S and b ^ T. In this situation we have 

l'{a) + l'{b) =l{a) -ai + l{b). 

Two cases arises: 

— l{a) + 1(b) — Weight{a, b) = ai. So 

l\a)+l'{b) = l{a)-ai+l{b) = l{a)-l{a)-l{b)+Weight{a,b)+l{b) = Weight{a,b). 

Hence, Ei C Ei/. 

— l{a) + l{b) — Weight(a,b) > ai. Obviously 

r{a) + r{b) = l{a) -ai + l{b) > Weight{a, b). 

□ 

Theorem 1. /// is feasible labeling and M is a Perfect matching in Ei, then 
M is a max-weight matching 0/. 

Proof. Suppose that M' is a perfect matching in G, since each vertex is 
incident to exactly one edge of M' we have: 

Weight{M')= ^ Weight{a,b) < J] l{v). 

{a,b)£M' t)6(AUB) 



So (AuB) K'^) is an upper bound for each perfect matching. Now assume 
that M is a perfect matching in Ef. 

Weight{M) = J2Ke) = 

eeM V&{AUB) 

It is obvious that M is an optimal matching. □ 
In the following, we briefly describe the basic Hungarian algorithm which 
computes the maximum many to many matching between two sets. The 
input bipartite graph G = {A U B,E) is a complete bipartite graph that in 
which \A\ = \B\ = n. 

Algorithm 1 The Basic Hungarian algorithm(A,i?) 
1: Initial > Find an initial feasible labeling / and a matching M in Ei 
2: Let l{bj) = 0, for all 1 < j < t 

3: Z(aj) = max*^]^ Weight{ai, bj) for all 1 < i < s 

4: M = 

5: while M is not perfect do 

6: Select a free vertex ai E A and set = {aj}, T = 
7: repeat 

8: while iVi(5) = T do 

9: Update the labels according to Lemma [1] 

10: Select 6, G A/'K^) - T 

11: if bj is not free then> {bj is matched to the vertex z, extend the 

alternating tree) 

12: S = SU z,T = TUbj. 

13: until bj is free 

14: Augment M 
return M 



In line 1, we label all points of B with zero and each point ai G A with 
max^^j^ Weight{ai, bj) to get an initial feasible labeling. Note that M can be 
empty. It is obvious that for computing the minimum cost many to many 
matching using the Hungarian algorithm we must weight the edge (oj, bj) by 
1/Weight{ai, bj). 



Lemma 2. Each augmenting path is a 4-vertex path. 

Proof. Suppose that the lemma is false. Let p = ai, bi, a2, &2, • • • , be an 
augmenting path with more than four vertices, that is k > 2. Note that ai 
and bk are free nodes. It is obvious that the first edge is in — M, so the 
second , third, and fourth edges of p are in M, E — M, and M, respectively. 
Since the third edge (02, 62) is in — M, the fourth edge (62, as) must be in 
M. Note that 62 is a free node. A contradiction. □ 

3. The algorithm 

In this section, we describe our new algorithm which is based on the well 
known Hungarian algorithm. Consider two point sets A = {oi, a2, . . . , a^} 
and B = {bi, 62, • • • , h} with s + t = n. Let Da = {ai, 0.2, ■ ■ ■ , as} and 
Db = /32, • • • , A} denote the demand sets of A and B, respectively. 
Let Ca = {o:'i,a'2, ■ ■ ■ ,a's} and Cb = {f^'i, (^'2^ ■ ■ ■ y P't} be the capacity 
sets of A and B, respectively. Without loss of generahty, we assume that 

Theorem 2. Let A and B be two sets with \A\ + \B\ = n, an MMDC match- 
ing between A and B can be computed in 0{n^) time. 

Proof. 

We first construct a bipartite graph as follows. Consider the complete 
bipartite graph G = {XUY,E) where X = AU A' and Y = B U B' U C (see 
Figured]). A complete connection between two sets is a connection that in 
which each element of one set is connected to all elements of the other set. 
We show each set of the vertices by a rectangle and the complete connection 
between them by a line connecting the two corresponding rectangles. 

Given A = {ai, 02, ... , a^} and B = {bi, 62, • • • , h}, there exists a com- 
plete connection between A and B such that the weight of (oj, bj) is equal to 
the cost of matching the point to bj for all 1 < i < s and 1 < j < t. Let 
B' = {b[, 62, ■ ■ ■ ,b[} and A' = {a[, a'g, . . . , a'^}, each point of A is connected 
to the all points of B' such that the weight of (a,, b'j) is equal to the weight 
of {ai,bj). There exists also a complete connection between the sets B and 
A' such that the weight of {a[, bj) is equal to the weight of (cj, bj). We have 
a set C = {ci, C2, . . . , Ch} that in which h = X]i=i (^i ~ X]j=i l^j- f^*^^' '^^ 
use C to get |X| = \Y\. Each vertex of A' is connected to all vertices of C 
with zero weighted edges. 
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Figure 1: Our constructed complete bipartite graph with h = '^i ~ X]'=i P'j- 

Now we apply our new algorithm, Algorithm [21 on above bipartite graph 
G. Let Cap{u) and Dem{u) denote the capacity and the demand of the 
vertex -u; so for all i, j we have Dem{ai) = ai, Dem{hj) = /3j, Cap{ai) = a[, 
and Capihj) = 

In our algorithm, a vertex x is free to another vertex y ii x is not matched 
with y in M and has at least one empty capacity. So Oj G A and a- G A' are 
called free vertices to a vertex b that are not matched with it in M, if 
Num{ai) < Dem{ai) and Num{a[) < Cap{ai) — Dem{ai), respectively. 
Also the vertices bj and b'j are free to another vertex that is not incident in 
M to them, when 

Num{bj) < Dem{bj) and Num{b'j) < Cap{bj) — Dem{bj), respectively. 

In fact, we save the current number of the vertices that are matched to 
the vertices of A, B, A', and B' in the arrays y4[l . . . s], B[l . . .t], A'[l . . . s], 
and B'[l . . .t], respectively; for example A[i] shows the number of the nodes 
that are matched to a^. The initial values of the arrays is 0; when a new 
point is matched to their representing node their values are increased by 1. 
Assume that Num{u) returns the number of the vertices that are matched to 
u so far. So Num{a.j) = A[i], Num{a^) = A'[i], Num{bj) = B[j], and finally 
Nu'm{b'j) = B'[j]. Note that the procedures IsFree{u) and IsMatched{u) 
return True if Num{u) < Cap{u) and Num{u) = Cap{u), respectively. So 
in the augmenting path a, b, c, d, a is free to b, b is matched to c, and d is 
free to c. Now we change the basic Hungarian algorithm as follows. 

We first label the vertices of our bipartite graph G using an initial feasible 
labeling in lines 2 — 4. Algorithm [2] has a while loop where O(n^) times 
iterates and ^l^i cti + ^]=i Pj edges are selected. In each iteration of our 
algorithm \M\ increases by 1. Let 

slcicky = min{/(a;) + l{y) - Weight{x,y)}. 



Algorithm 2 The MMDC Hungarian algorithm(DA, CA, DB, CB) 
1: Initialize > Find an initial feasible labeling / and a matching M in Ei 
2: Let l{hj), l{h'-) = 0, for alll<j <t 

3: l{cLi) = m^Bx^j^iijiiayiiW eight{ai, bj), Weight{ai, bj)) for all 1 < i < 

s 

4: l{a'j) = max*^^ Weight{a[, bj) for all 1 < i < s 

5: Let M = 

6: while {u e AU A', with IsFree{u)} 7^ do 
7: Select u e AU A' with IsFree{u) 
8: Set^ = {M},T = 

9: repeat 

10: while Ni{S) = T do 

> Update the labels according to Lemma [1] 
11: Let ai = mms^(,s,tj<^T{l{si) + l{tj) - Weight{si, tj)} 

{l{v) — ai if V E S 
llv) + ai if V eT 
l{v) Otherwise 

13: Select y e Ni{S) - T 

14: if IsMatched{y) then > {Num{y) = Cap{y)) 

> {y is matched to some vertices z) 
15: S = SU{z\{z,y)eM},T = TU{y}. 

16: until IsFree{y) 



17: 



Augment{M) 



In line 17 ofAlgorithm m the values of all slacks must be updated when a 
vertex is moved form S to S. This is done in 0{n) time. During our algorithm 
s + t = n vertices are moved from 5^ to S, so it takes the total time of Oln"^). 
In lines 11, we can compute the value of ai by: 

ai = min slacky, 

in 0{n) time. After computing the value of ai and updating the labels of 
the vertices, we must also update the values of the slacks. This can be done 
using: 

^ Tslacky = slacky — ai. 

In each iteration the value of ai may be computed at most 0{n) times, that 
takes 0{n) time each time, so running each iteration takes at most O(n^) 
time. Our algorithm has 0{n^) iteration with 0{n'^) time, so it runs in O(n^) 
time. 

□ 

4. Conclusion 

In this paper, we presented an O(n^) time and 0{n) space algorithm for 
computing an MMDC matching between A and B with total cardinality n. 
In fact, we modified the basic Hungarian algorithm to get a new algorithm, 
called the MMDC matching algorithm. Then, we construct a bipartite graph 
G and apply our new algorithm on G. 
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